2011 Fermi Symposium, Roma, May 9-12 



First Detection of a Pulsar above 100 GeV 

A. N. Ottc 

Santa Cruz Institute for Particle Physics (SCIPP), 
Physics Department Univsersity of California Santa Cruz, 
1156 High Street, Santa Cruz, CA 95064, USA 

for the VERITAS Collaboration 
http ://veritas.sao. arizona. edu / 

We present the detection of pulsed gamma-ray emission from the Crab Pulsar above 100 GeV 
with the VERITAS array of atmospheric Cherenkov telescopes [l| . Gamma-ray emission at theses 
energies is not expected in present pulsar models. We find that the photon spectrum of pulsed 
emission between 100 MeV and 400 GeV can be described by a broken power law, and that it 
is statistically preferred over a power law with an exponential cut-off. In the VERITAS energy 
range the spectrum can be described with a simple power law with a spectral index of -3.8 and a 
flux normalization at f 50 GeV that is equivalent to 1% of the Crab Nebula gamma-ray flux. The 
detection of pulsed emission above 100 GeV and the absence of an exponential cutoff rules out 
curvature radiation as the primary gamma-ray-producing mechanism. The pulse profile exhibits 
the characteristic two pulses of the Crab Pulsar at phases 0.0 and 0.4, albeit 2-3 times narrower 
than below 10 GeV. The narrowing can be interpreted as a tapered particle acceleration region in 
the magnetosphere. Our findings require that the emission region of the observed gamma rays be 
beyond 10 stellar radii from the neutron star. 



I. INTRODUCTION 

One of the most powerful pulsars in gamma rays is 
the Crab Pulsar %% PSR J0534+220, which is the 
remnant of a historical supernova that was observed 
in 1054 A.D. It is located at a distance of 6500 light 
years, has a rotation period of «33 ms, a spin-down 
power of 4.6 x 10 38 erg s _1 and a surface magnetic field 
of 3.78 x 10 12 G 0. 

Within the corotating magnetosphere, charged par- 
ticles are accelerated to relativistic energies and 
emit non-thermal radiation from radio waves through 
gamma rays. In general, gamma-ray pulsars exhibit 
a break in the spectrum between a few hundred MeV 
and a few GeV. Mapping the cut-off can help to con- 
strain the geometry of the acceleration region, the 
gamma-ray radiation mechanisms and the attenuation 
of gamma-rays. Based on previous measurements and 
present theoretical understanding that the dominant 
gamma-ray emission mechanism is curvature radia- 
tion, it is widely believed that the shape of the spectral 
break is best described by an exponential cut-off. 

Although measurements of the Crab Pulsar spec- 
trum are consistent with a power law with exponen- 
tial cut-off, flux measurements above 10 GeV are sys- 
tematically above the best-fit model, suggesting that 
the spectrum is indeed harder than a power law with 
exponential cut-off However, the statistical un- 

certainty of the previous data was insufficient to allow 
a definite conclusion about the spectral shape. In this 
paper we summarise the recent detection of the Crab 
Pulsar above 100 GeV with VERITAS that rules out 
that the spectrum above the break is described by an 
exponential cut-off. In Section 2 we describe the ob- 



servation and analysis. The results are presented in 
Section 3 and we close the paper with a discussion in 
Section 4. 

II. OBSERVATION AND ANALYSIS 

VERITAS, the Very Energetic Radiation Imaging 
Telescope Array System, is an array of four 12 m di- 
ameter imaging atmospheric Cherenkov telescopes lo- 
cated in southern Arizona, USA After evidence 
for pulsed emission was seen in 45 hours of data from 
the Crab Pulsar that were recorded between 2007 and 

2010, a deep 62-hour observation was carried out on 
the Crab Pulsar between September 2010 and March 

2011. The observations were made in wobble mode 
with a 0.5 degree offset. After eliminating data taken 
under variable or poor sky conditions or affected by 
technical problems, the total analysed data set com- 
prises 107 hours of observations (97 hours dead-time 
corrected) carried out with all four telescopes. The 
data were taken with the standard VERITAS trig- 
ger setting, and analysed with the standard VERITAS 
analysis tools. 

A. Event Reconstruction 

In the analysis, the recorded atmospheric shower 
images are processed with a standard moment anal- 
ysis Q and the energy and arrival direction of the 
primary particle are calculated Q. After event re- 
construction, a selection is performed to reject events 
caused by charged cosmic rays. The selection crite- 
ria were optimised a priori for highest sensitivity by 
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FIG. 1: VERITAS pulse profile of the Crab Pulsar at > 120 GeV. The shaded histograms show the VERITAS data. The 
pulse profile is shown twice for clarity. The dashed horizontal line shows the background level estimated from data in the 
phase region between 0.43 and 0.94. The data above 100 MeV from the Fermi-LAT @] are shown beneath the VERITAS 
profile. The vertical dashed lines in the panels mark the best-fit peak positions of PI and P2 in the VERITAS data. 
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FIG. 2: Enlarged view of the main pulse or PI of the Crab 
Pulsar pulse profile. See text. 



FIG. 3: Enlarged view of the interpulse or P2 of the Crab 
Pulsar pulse profile. See text. 



assuming a simple power-law energy spectrum for the 
Crab Pulsar with an index a = —4 and a flux nor- 
malisation of a few percent of the Crab Nebula flux 
at 100 GeV. The analysis threshold is 120 GeV. 

For the pulsar analysis, the arrival times of the se- 
lected events are transformed to the barycenter of the 
solar system. After barycentering, the phase of the 



Crab Pulsar is calculated for each event using con- 
temporaneous ephemerides of the Crab Pulsar that 
are published monthly by the Jodrell Bank telescope 
i- 

The results were confirmed by a separate analysis of 
the data made using an independent analysis package. 
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III. RESULTS 
A. Pulse Profile 

The phase-folded event distribution, hereafter pulse 
profile, of the selected VERITAS events is shown in 
Figure [U The most significant structures are two 
pulses with peak amplitudes at phase 0.0 and phase 
0.4. These coincide with the locations of the main 
pulse and interpulse, hereafter PI and P2, which are 
the two main features in the pulse profile of the Crab 
Pulsar throughout the electromagnetic spectrum. In 
order to assess the significance of the pulsed emission, 
we use the H-Test [1(|. The test result is 50, which 
translates into a statistical significance of 6.0 standard 
deviations that pulsed emission is present in the data. 

The pulse profile has been characterised by an un- 
binned maximum-likelihood fit; see the solid black line 
in Figures [5] and El In the fit, the pulses are modeled 
with Gaussian functions, and the background is de- 
termined from the events that fall between phases 0.5 
and 0.9 in the pulse profile (referred to as the off-pulse 
region) . 

The positions of PI and P2 in the VERITAS 
data are thus determined to lie at the phase values 
-0.0026 ± 0.0028 and 0.3978 ± 0.0020, respectively 
and are shown by the vertical lines Figure 1. The 
full widths at half maximum (FWHM) of the fitted 
pulses are 0.0122 ± 0.0035 and 0.0267 ± 0.0052, re- 
spectively. The pulses are narrower than those mea- 
sured by Fermi-LAT at 100 MeV by a factor of two to 
three. The energy- dependent narrowing of the pulses 
is a strong probe of the nature of the magnetospheric 
particle acceleration region and can be used to shed 
some light on its geometry, electric field, and gamma- 
ray emission properties. If gamma rays observed at 
the same phase are emitted by particles that prop- 
agate along the same magnetic field line [Tljj , then a 
possible explanation of the observed narrowing is that 
the region where acceleration occurs tapers towards 
the neutron star. However, detailed calculations are 
necessary to explain fully the observed pulse profile. 



0.421 for P2. This is the ±2 standard deviation inter- 
val of each pulse as determined from the maximum- 
likelihood fit. The significance of the excess in these 
regions is 4.7 standard deviations for PI and 7.9 stan- 
dard deviations for P2. The ratio of the excess events 
and thus the intensity ratio of P2/P1 is 2.4 ± 0.6. If 
one assumes that the differential energy spectra of PI 
and P2 above 25GeV can each be described with a 
power law, dN/dE cx E a , and that the intensity ratio 
is exactly unity at 25 GeV [j| , then the spectral index 
a of PI must be smaller than the spectral index of P2 
by ap2 — ap^ = 0.56 ± 0.16. 



C. Phase-averaged spectrum 

The gamma-ray spectrum above 100 GeV was mea- 
sured by combining the signal regions around PI and 
P2 defined above. This can be considered a good 
approximation of the phase-averaged spectrum since 
no "bridge emission", which is observed at lower en- 
ergies, is seen between PI and P2 in the VERI- 
TAS data. However, the existence of a flux com- 
ponent that originates in the magnetosphere and is 
uniformly distributed in phase cannot be excluded 
and would be indistinguishable from the gamma-ray 
flux from the nebula. Figure 0] shows the VERI- 
TAS phase-averaged spectrum together with measure- 
ments made with Fermi-LAT and MAGIC. In the en- 
ergy range between 100 GeV and 400 GeV measured 
by VERITAS, the energy spectrum is well described 
by a power law dN/dE = Ax (E/150 GeV)", with 
A = (4.2 ± 0.6 stat + 2.4 syst - 1.4 syst ) x KT 11 TeV" 1 
cm~ 2 s _1 and a = — 3.8±0.5 s ^ a ^±0.2 S y S ^. The detec- 
tion of pulsed gamma-ray emission between 200 GeV 
and 400 GeV, the highest energy flux point, is only 
possible if the emission region is at least 10 stellar 
radii from the star's surface 113. 



D. The spectral energy distribution between 
100 MeV and 300 GeV 



B. Spectral differences between PI and P2 

Along with the observed differences in the pulse 
width, the amplitude of P2 is larger than PI in the 
profile measured with VERITAS, in contrast to what 
is observed at lower gamma-ray energies where PI 
dominates (see Figure [1]). It is known that the ra- 
tio of the pulse amplitudes changes as a function of 
energy above 1 GeV [|[ and becomes near unity for 
the pulse profile integrated above 25 GeV @ . 

In order to quantify the relative intensity of the two 
peaks above 120 GeV, we integrate the excess between 
phase -0.013 and 0.009 for PI and between 0.375 and 



Combining the VERITAS data with the Fermi- 
LAT data we can place a stringent constraint on 
the shape of the spectral turnover. The previously 
favoured spectral shape of the Crab Pulsar above 
1 GeV was an exponential cut-off dN/dE = A x 
(E/E ) a exp(—E/Ec), which is a good parametrisa- 
tion of the Fermi-LAT Q and MAGIC @ data. We 
note that the Fermi-LAT and MAGIC data can be 
equally well parametrised by a broken power law but 
those data are not sufficient to distinguish significantly 
between a broken power law and an exponential cut- 
off. The VERITAS data, on the other hand, clearly 
favour a broken power law as a parametrisation of 
the spectral shape. The fit of the VERITAS and 
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FIG. 4: Spectral energy distribution (SED) of the Crab Pulsar in gamma rays. VERITAS flux measurements are shown 
by the solid red circles, Fermi-LAT data [j| by green squares, and the MAGIC flux point [f| by the solid triangle. The 
empty symbols are upper limits from CELESTE, MAGIC, and Whipple. The bow tie and the enclosed dotted line give 
the statistical uncertainties and the best-fit power-law spectrum for the VERITAS data using a forward-folding method. 
The result of a fit of the VERITAS and Fermi-LAT data with a broken power law is given by the solid line and the 
result of a fit with a power-law spectrum multiplied with an exponential cut-off is given by the dashed line. Below the 
SED we plot x 2 values to visualise the deviations of the best-fit parametrisation from the Fermi-LAT and VERITAS flux 
measurements. 



Fermi-LAT data with a broken power law of the form 
A x (E/E ) a /[1 + (E/E ) a ~P] results in a X 2 value 
of 13.5 for 15 degrees of freedom with the fit param- 
eters A = (1.45 ± 0.15 stat ) x 1(T 5 TeV" 1 cm~ 2 s _1 , 
E = 4.0 ± 0.5 stat GeV, a = -1.96 ± 0.02 stat and 
j3 = — 3.52±0.04 8 ^ a ^ (see solid black line in Figure 2). 
A corresponding fit with a power law and an expo- 
nential cut-off yields a \ 2 value of 66.8 for 16 degrees 
of freedom. The fit probability of 3.6 x 10 -8 derived 
from the \ 2 value excludes the exponential cut-off as 
a viable parametrisation of the Crab Pulsar spectrum. 

IV. DISCUSSION 

The detection of pulsed gamma-ray emission above 
100 GeV provides strong constraints on the gamma- 
ray radiation mechanisms and the location of the ac- 
celeration regions. For example, the shape of the spec- 
trum above the break can not be attributed to curva- 
ture radiation because that would require an expo- 
nentially shaped cut-off. Assuming a balance between 
acceleration gains and radiative losses by curvature 
radiation, the break in the gamma-ray spectrum is 
expected to be at Eb r = 24 GeV ?7 3 ^ 4 V?, where r\ is 
the acceleration efficiency {r\ < 1) and £ is the ra- 



dius of curvature in units of the light-cylinder radius 
[i~3l ]. Though £ can be larger than one, only with an 
extremely large radius of curvature would it be possi- 
ble to produce gamma-ray emission above 100 GeV 
with curvature radiation. It is, therefore, unlikely 
that curvature radiation is the dominant production 
mechanism of the observed gamma-ray emission above 
100 GeV. Two possible interpretations are that either 
the entire gamma-ray production is dominated by one 
emission mechanism different from curvature radia- 
tion or that a second mechanism becomes dominant 
above the spectral break energy. 
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